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Abstract — This paper investigates decentralized control for 
linear time-invariant systems with norm-bounded time-varying 
structured uncertainty. We derive LMI conditions ensuring 
robust stability. We also show how to incorporate performance 
robustness, where performance is measured by the H„ gain of 
the decentralized system. The potential of the proposed 
technique has been demonstrated by a simulation example of 
three coupled inverted pendulums. 

Index Terms —Decentralized control, LMI, norm-bounded 
uncertainty, robust stability. 


I. Introduction 

A large scale system is often decomposed into a set of 
inter-connected subsystems to enable control implementation. 
Thus, research on decentralized control has been attracted the 
attention of many researchers since the seventies. A sufficient 
condition [1] for connective stability in terms of bounds of 
uncertain inter-connections between subsystems was studied. 
This connective stability laid the foundation of decentralized 
control. In [2] and [9]-[10] considered LTI systems with 
time-varying but bounded parametric uncertainties with the 
assumption that the uncertainties are in the range space of 
input matrix (matching condition). In [3] and [8], the 
matching conditions of the system are not valid. Thus the 
authors decomposed the uncertainties into matched and 
mismatched parts. The matched portion can be dealt with by 
feedback control while the mismatched portion was 
constrained by norm conditions in terms of state signals. In [4], 
the system with time-varying, norm-bounded parametric 
uncertainties with no matching assumption were studied. 
Sufficient conditions for stability were established in terms of 
a scaled H x control system and a robust II XJ decentralized 
controller was designed accordingly. 

Note that decentralized control design is, in general, more 
difficult than that of centralized control. The difficult arises 
due to the fact that the controller for each subsystem is 
restricted to only using the local state or measurements in 
order to stabilize (with prescribed performance) the overall 
system. As in most control systems, a mathematical model 
usually cannot describe a dynamic system exactly. Also, a 
dynamical system is mostly working in a changing 
environment. This is especially true for large scale systems. 

The uncertainties considered in this paper consist of real 
time-varying norm-bounded structured uncertainties that are 
realized by a linear fraction representation (LFR) in state 
space framework. LFR has gained significant attention due to 
its structured representation in terms of matrices [5]. The LFR 


consists of an LTI system, connected with a diagonal 
feedback element depending on the setting of uncertainties. 
No “matching condition” is imposed on the uncertainties 
because of the system representation. We present LMI-based 
system analysis and controller synthesis in which Lyapunov 
matrix inequalities in terms of constraints can be cast as a 
convex optimization problem [7], [8]. 

The paper’s outline is as follows. In Section 2, we 
introduce the notation used in the paper and describe the 
system to be analyzed and controlled. Also, we define the 
decentralized control problem considered in the paper, 
including the stability and performance robustness. The 
performance is measured by L 2 gain. Section 3 is devoted to 
the analysis of the open loop system, based on the robustness 
requirement. The tools for analysis are Lyapunov stability 
theory and convex optimization techniques. In Section 4, we 
address the corresponding controller synthesis problems for 
static state-feedback. In Section 5, we illustrate our methods 
on three inverted pendulums control problems. 


II. Problem Formulation 


A. Notation 


Throughout the paper the subscript “i” denotes zth 
subsystem. I denotes the identity matrix with its size 
determined from context, q denotes the rth element of the 
vector q and diagM./l) is defined as 


diag (A, B ) = 


A 

0 


0 

B 


We reserve the subscript “d” to express 


A d := diag (A,), =1 


A, e R 


tij xrij 


The subscript “D” is for 

B d :=diag (S ; ) i=1 „ B, e 


B. Problem to Solve 

The purpose of this subsection is to define a framework on 
which our approach to decentralized control is based. 
Consider a class of linear time-invariant continuous-time 
systems E with norm-bounded time-varying structured 
uncertainty, composed of n subsystems Z,. 

n 

X:= U 'Li(x,u i , p t ,Wi) 

i =1 

where x is the state of system X and is defined as 

[ t t tY D n , 

x = [ x i ■■■ x i x„ \ . x t e R 

The subsystem Z,- is defined as follows. 
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x i = A i x i + I A i x j + B u,i u i + BpjPi + B w,i w h 
Zi — C^jXj + D zu jUj + D~^i Wj , 


9i = CgjXf + x C J qi Xj + D qu j Ui 
j=hj*i 


(1) 


" DqpjPi Dqw,i w i' 


Pi = (t)qj > s[(t) 


<1 


r = 1, 2 ,...., 11 , 




where x, denotes the state of subsystem i , xj is the state that 
acts on subsystem i from subsystem j, u, is the decentralized 
control input, signal w, is the external inputs which include 
any unknown disturbance to be rejected as well as the 
reference command to be tracked. z„ the controlled variable, 
may include tracking error or a cost of the input m,. Assuming 
that the system is full state-feedback, that is y ; = x t . The 
matrices A, B, C, and D with all superscripts and subscripts 
shown in (1) are assumed to be real constant matrices with 
appropriate dimensions. If x, has no influence on x„ the 
associated system matrix Af will be set to zero. Notice that /?,■ 
and q , represent the output and input vectors of system 
perturbations while the pf and q[ correspond to the rth 
element of p, and q, respectively. 

J „r » 


Pi = 


: R 


9i = 


1 

9i 


r 

9i 



q- e R 


5/(0 is the real valued unknown but bounded uncertainty and 
it is assumed to be a time-varying scalar quantity. Let 

Ai d ■= dia g(#'(0) r=1; .. . in? ., $i (0 e R. 

It is assumed that the uncertainty A,- d belongs to a compact set 

J&A := K'.rf : |[ A i,d |j - 1’ 

where the prefix stands for “ball” or “bound”. Notice that 
the above-mentioned bounds applied in (1) can also be 
replaced by pj p, < qj q,. The state feedback closed-loop 
system of (1) is shown in Fig. 1 which mainly gives one the 
idea and meaning of how the system input, output, and 
uncertainty are represented. In Fig. 1 the labeled Z, represents 
the open-loop interconnection and contains all of the known 
elements including the nominal plant model, performance and 
the uncertainty weighting functions. The A l rJ block contains 
all the uncertain elements and is assumed to be within a 
compact set jg A which parametrizes all of the parametric 
uncertainty in the problem. There are four sets of inputs 
entering into Z,: perturbation output /?„ disturbance w„ time 
varying state inputs from other subsystem x, and /V i, and 
decentralized controls «,(f). Three sets of outputs are: 
perturbation input q h controlled variable z r system state x„ 
and state feedback gain K,. 

The robustness specifications we address in this paper are 
that the L 2 gain of input/output signals for each individual 
subsystem is less than a specific number y,, which corresponds 
to the peak gain of the freq. response G,(joJ), i.e. 

|G i (./®)|| oo =supff(G i ( i /®)) 

CO 


Structured Uncertainty 



Fig. 1 Closed-loop for subsystem i 


where IIG/I*, represents the II X norm of transfer function 
matrix from disturbance w, to controlled variable z,. Thus, 
according to the system described in (3.1) and the required 
robustness properties the following robust decentralized 
control problem is addressed. 

Find a decentralized feedback controller w,(t) as a linear 
function of local state vector x,(t), i.e. u,{ t) = K,x,(X) in which Kj 
is the feedback gain for subsystem i such that the resulting 
closed-loop system is robustly stable with respect to 
decentralized //„ disturbance attenuation for all admissible 
perturbations which lie in a compact set B v Sufficient 
conditions are 

(1) Closed-loop system is asymptotic stable. 

(2) IIG,IL<y,, 


III. System Analysis 

One of the powerful methods to analyze system stability is 
Lyapunov stability theory. Although there is no trivial method 
to establish a Lyapunov function, by experience, among 
Lyapunov functions that have been suggested and extensively 
used, the quadratic Lyapunov functions, i.e. V(£) = c'Lz, have 
been proved to be efficient and easily implemented [7]. 

Let’s consider an equivalent system of (1) 

n 

x i = A i x i + Z A i' x j + BpjPi + B w,i w i , 
i=hj*i 

— C~jXj + D^yyjWj, ( 2 ) 

9i = Cq,i x i DqpjPi Dqw,i w i’ 

T s T 
Pi Pi < 9i 9i- 

which is denoted by IT/. Therefore, we have the Theme 1. 


Theme 1 The overall system n := U” = iT~L(jq) is robust 


decentralized stable with respect to H x disturbance 
attenuation if there exists a symmetric matrix L d such that the 
following matrix inequalities are satisfied 
L d > 0, A rf >0, F rf >0, and 


BLK11 BLK12 BLK13 
(BLK12) 7 BLK22 0 
(BLK13) T 0 BLK33 


<0 


(3) 


where 
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BLK11 = 


( A d L d + L d A d + L d B a, D B a,D L d ^ 
+ Cz.D C z,D +C q A d C q + ^d 


BLK12 = L d B p D + Cq A d D qp D 
BLK22 = Dg p D A d D qp D -A d 
BLK13 = L d B w D + C ! z 1)0 zw p 

BLK33 = Dl wD D zwI) -T d . 

Proof. See appendix A. 

Remarks: 

1) The last inequality in (3) is a quadratic inequality in matrix 
variable L d can be further decomposed into LMI by Schur 
complement [7], [8] as shown in (4). 

BL11 BLK12 BLK13 L d B a D 

\T 


(BLK 12) BLK22 0 

T 0 BLK33 

0 


(BLK13) 

_( L d B a,D) T 
where 


0 


0 

0 

-I 


<0 


(4) 


-)• 


Qd 


(QBLK13) 

where 
QBLK11 = 


0 

0 


QBLK33 0 
0 -H 


f Qd A d +A dQd +B a,D B a,D 


<0 (5) 


+ B p,D M d B ]>,D + B w,D N d B w,D 


QBLK12 = B p D M d D T qP ' D + (Q d c' q , 
QBLK 13 =B wD N d Dl wD +Q d C^ D , 


QBLK22 =-\M d ~D qp D M d D^ D 


BLK22 0 
0 BLK33 


<0 


or, equivalently, 

QBLK22 0 
0 QBLK33 


<0 


( 6 ) 


Maximize Trace N d 

Subject to Q d = Qj > 0, N d > 0, M d > 0, and (5) with 
well-posed condition (6). 


IV. Controller Synthesis 

Consider the subsystem X, with static state-feedback 
controller M,(t) = KjXj( t). The decentralized closed-loop 
system can be written as 


' X A i •*j + B p d Pi + B w,i w i 


(7) 


BL11 - \A d L d +L d A d + C zD C zD + C q A d C q + Q, 

2) A dual condition of (4) can be established by a change of 
variable. Let Q d - L r / , the condition (4) can be 
equivalently written as 

QBLK11 QBLK12 QBLK13 Q d 


x i ={A\ +B ui K i )x i + X aJ 
j=hj*i 

+ D Z ujKi)xi + 

<li = ( C qJ + D qu,i K i ) x i + I C q,i x j 

DqpjPi "*■ Dqw,i w i 

Thus, the overall closed-loop system may be written as 
x = ( A d + B u,D K D ) x + B a.D x + B p,D P + B w,D w 
z = ( C z,D + D zm.D K d) x +D zy\> ,D W ( 8 ) 

9 = ( C q + D qu,D K D) x+ D qp.DP + D 'qw.D w 
In the view of Theme 1 and its remarks, the sufficient 
condition for system (7) to have robust decentralized H„ 
disturbance attenuation is 
Q d > 0, M d > 0 , N d > 0, and 

KBLK11 KBLK12 KBLK13 


(KBLK12 y QBLK22 


(KBLK13) 

Qd 


o 

QBLK33 

0 


Qd 

0 

0 

-H d 


<0 


(9) 


where 


KBLK11 = 


Qd( A d +B u.D K d) +( A d + B H,D K o)Qd 

+ B a,D B a,D + B p,D M d B ]i,D + B w,D N d B w.D 


QBLK33 = -\I-D ZWiD N d D 1 WtD ), 

Al d = A/ 1 , N d = IV' , and H d = Q/ 1 . In fact, decentralized 
stability with performance requirement can be implied by 
the existence of Q d > 0 satisfying N d > 0, M d > 0, and (5). 
3) The well-posed condition for (2) is incorporated in (4) or 
(5), since it implies 


KBLK12 - B pD M d D qpD + Q d (C q + D qu D K D ) T 
KBLK13 = B wD N d D^ wD +Q d (C zD + D zll D K D ) T . 


Note that the last inequality in (9) is not convex since it jointly 
depends on Q d and K D . By the change of variables Y = K D Q d , we 
obtain 


YBLK11 
(YBLK12 j 1 
(YBLK13) 1 


Qd 


where 


YBLK11 = 


YBLK12 YBLK13 Q d 
QBLK22 0 0 

0 QBLK33 0 
0 0 -H d 

,T D T 


<0 ( 10 ) 


qu,D 


4) From Theme 1 and Remark 1) and 2), we see that a 
sufficient condition to determine system stability and 
performance can, in turn, be cast into convex optimization 
problem that is shown in [7], [8], The problem can be 
easily formalized as an LMI feasible problem or can be cast 
into eigenvalue problem (EVP) as follows. 


Qd A d +Y B u.D + A dQd +B u.D Y 

+ B a,D B a.D + B p,D M d B ]i,D + B w,D N d B w,D 
YBLK12 =( B pD M d D T qpD +Q d C T q +Y T D i 

YBLK13 = {b k/) N d D zw j ) + Q d cl D +V T 

(10) is LMI in matrix variables Q d , Y, M d , and N d and thus can 
be cast into a convex optimization problem as follows. 
[Maximize Trace N d 


Subject to Q d > 0, N d > 0, M d > 0, and (10). 


(ID 
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Once the maximized solutions are found for the problem, then 
a stabilizing state feedback gain is obtained by 

K D =YQj l (12) 

V. Numerical Examples 

Consider the system in Fig. 2 consisting of three inverted 
pendulums of point masses and length /,. The pendulums 
interact via three springs and three dampers of stiffness k, t and 
damping i, j = 1,2,3, and fej. The distances from attached 
point of springs and dampers to the platform baseline are a,. 



Fig. 2 Three Inverted Pendulums Systems 

In this example, we will demonstrate step response for the 
system using LMI-based decentralized control and show the 
flexibility of LMI-based approach. In order to achieve zero 
steady state error, i.e. e(t^oo) = 0, an integral control will be 
added as an extra state for each subsystem. The system 
dynamics for pendulums are written in general form. 


0i= — 
m.; 


3 

X 

j=hj*i 


bij&i + —r k ij & i 

h 


l, 


/ b -A + V k u°j 

V 1 l i J 

+ BL 0i+ Lu t 

h l; 


\ 


(13) 


Multiplicative uncertainties are used as shown below. 


nil — w,(l+A m ,5 mi ), bjj — h_jj{ 1 +A/ ) ^(5/ p; y), k,j — k^i 1 +A^5^) 


in which A with subscripts, as uncertainty weighting, is the 
percentage perturbation from nominal value of each 
parameter and 5 with subscripts varies from -1 to 1, i.e. -1 < 5 
< 1. Plant data and weighting of uncertainties are shown in 
Table 1. Thus, the system can be written as (8). In this 
example, we let D zuD =0 to demonstrate a singular case in 
traditional //„ setup while the control saturation can still be 
enforced as constraints shown in [7]. A full state-feedback 
static gain can be obtained by computing feasible solution of 
(11). If feasible solution does exist, then static gain can be 
obtained by (12). The simulation was done using Matlab with 
LMI Control Toolbox and static gain 


K d = 


K t 
0 


0 

k 2 


0 0 K 


where 0 represents the zero vector of appropriate dimension. 


Table 1 



rrij 

m 2 

m 3 

bn 

bn 

b 2 3 

nominal 

i 

1 

i 

1 

.51 

1 

weighting 

.14 

.29 

.19 

1 

1 

1 


kn 

ki3 

k 2 3 

,?=io 

nominal 

1 

.51 

1 

h= 

, 6=1.2 

(3=1.1 

weighting 

1 

1 

1 

Q, i —— &3 

=0.5 


Ky = [-25.0328 -8.0267 -7.2189] 

K 2 =[—28.4028 -10.0355 -7.5556], and 
/C 3 = [-25.9417 -8.4924 -7.3352]. 

We obtain the overall system poles at 

-0.9462 + y'0.2215 , -7.6443, -1.0111, -8.6303, -0.7215, 

-7.3363, and -0.9470+jO. 1099. 

Notice that the uncertain data set in Table 1 can be used to 
form a data convex hull. The extreme points (or vertices) of 
the convex hull are taken from the corresponding 5 values, i.e. 
5=1 and 8=-l. We assume the system uncertainties vary 
within this convex hull. It is obvious that the mass weighting 
in the Table 1 should not be equal to 1 to avoid the possibility 
that zero mass occurred in the denominator of (7). From the 
chosen data in the Table 1 we demonstrate the simulation 
results by varying the uncertainties from one vertex, which the 
stiffness and damping are twice of their nominal values, to 
another vertex, which the corresponding stiffness and 
damping are zero values. This can be done by switching 5 
from 1 to -1. The sudden humps at time 5 sec in Fig. 3 and Fig. 
4 show the results due to this change. Notice that it is an 
instantaneous switch. Next, we turn our attenuation to 
examine conservativeness of our design. This can be checked 
by expanding the vertex of convex hull to see how large it can 
be to cause instability while keeping the same control gain. To 
examine every vertex of the convex hull is tedious, thus we 
will only expand the vertices shown in the simulation as 
before to get idea of conservativeness. The simulation results 
show the weighting of stiffness and damping can be 
simultaneously increased to approximately 3 which shows the 
design is conservative. In fact, the results are not surprised 
because Lyapunov function provides only sufficient condition 
to stabilize a system, which in turn causes a conservative 
design. 
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Fig. 4 Control inputs 


Appendix 


To prove Theme 1, we will first introduce the following 
identity and define the notation that will be used in the proof. 
Identity: 

A T B T CD + D T C T BA 

= —{D - C T BA) T (D - C T BA) + D t D + A T B T CC T BA. 
Notations: 

The diag(») and the usage of subscripts “<£” and “D” have 
been defined in Section II. The following notations will be 
used throughout the paper. 

A,- = diag (^|\ = \ ... nqi , A d = diag(A / ) ;=1 ... ;) , 

: dia g k? 1=1,- -n ’ C z.D = diag). =1 n , 


r, =< 


a i = £l . =diag(a,). = | _ n , 

k=l,k*i 


B. 


/; =diag(/4,) (=| n , and 

,=k ■■■ A -1 0 a' +1 A"J- 


B a ,i = [A 

Proof of Theme 1 : 

We consider the quadratic Lyapunov function V,(x,)= 
x, r L,x,. The consideration of performance and perturbation 
will be incorporated in the derivative of Lyapunov function 
for stabilization of overall system. We have 


z 

1=1 


VfXjf + zJZi-rfwjwi 


n qi T T 

, -dr qr, r r r r\ 

+ Y.A i (q i qj-Pi Pi) 


<0 


(14) 


V r =1 

Substituting (2) into (14), applying the identity, and 
rearranging the expression, we obtain the inequality as 
follows, 

f t’ A 


n 

z 

(=1 


Xj Zl\Xi + + 






T 


r 1 

T 



Xj 


x i 

X 


X 



• 



S 2 


+ 

Pi 

S 3 

Pi 

_Pi _ 


Pi 








_Wi_ 


_Wj_ 


■T 


Z ( Xj-A / LjXiYixj-A / A'X/) 

V f=l. i*i j 


<0, (15) 


where 


jT 

= A M' t ' "*■ C 


^ h ■ -T^ 

LA/A/ 

U=y 1 


Lj + ocfl + C' z i C‘ z j, 


■—< Z 
-2 


C T A C ■ C T A D 


^ qj ± *-i ± ^ qp,i 

T T 

qpd L i Lq i D qp jAiD qp i — 

0 LiB i @i 


and 


«;,,a o o 


o ©', 


(eif 


for 0j - Lj-Sw.j + ^nv,/ 


and & 2 =Dl v i D zw i - r /l . 


By summing up (15), we obtain the inequality 


X 


BLK11 

BLK12 

BLK13" 

X 

p 


(BLK12 ) T 

BLK22 

0 

p 

w 

T 

(BLK13) T 

0 

BLK33 

w 


n n -T -j, ;T 

-z Z (Xj-A/ LiXifixj-A/ L iXi ) 

v i=ij=hi*i 

A sufficient condition for the above inequality to hold 
requires Theme 1. This completes the proof. 
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